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ABSTRACT 

The observational and theoretical status of the search for missing cosmological baryons is sum- 
marized, with a discussion of some indirect methods of detection. The thermal interpretation of the 
cluster soft X-ray and EUV excess phenomenon is examined in the context of emission filaments, which 
are the higher density part of the warm hot intergalactic medium (WHIM) residing at the outskirt 
of clusters. We derived an analytic radial profile of the soft excess surface brightness using a simple 
filament model, which provided us a means of observationally constraining the WHIM parameters, 
especially the total mass budget of warm gas associated with a cluster. We then pointed out a new 
scenario for soft excess emission, viz. a cluster that can strongly lens the soft X-rays from background 
WHIM knots. If, as seems quite likely, the missing baryons are mostly in the WHIM halos of galaxy 
groups, the lensing probability will be quite high (~ 10 %). This way of accounting for at least part of 
a cluster's soft excess may also explain the absence of O VII absorption at the redshift of the cluster. 
Subject headings: galaxies: clusters; cosmology: large-scale structure of universe 



1. INTRODUCTION 

The location of the 'missing baryons' in the Universe 
is an open and important question of contemporary cos- 
mology, as interesting as it is ironic, because the prob- 
lem manifests itself as a deficit in the mass budget which 
arises only at low redshifts, i.e. in the space near us. Ob- 
servationally the total baryonic content in stars, galaxies, 
and clusters of galaxies {fib = {2.1^^'^)hQ^ %, Fukugita 
et al 1998) is only about half of the amount required by 
Big Bang Nucleosynthesis models (fif, — (3.9 ± 0.5)hQj 
%, Buries & Tytler 1998) or from measurements of the 
cosmic microwave background {fib — (4.6 ± 0.2)/i(^7, Ko- 
matsu et al. 2008 WMAP5, consistent with the Ben- 
nett et al 2003 WMAPl and Spergel et al 2006 WMAP3 
results). Cosmological hydrodynamic simulations have 
shown that this missing 50 % of baryons is concealed in 
a tenuous filamentary gas of temperature 10^ - 10^ K, 
currently referred by many to as the WHIM (the Warm 
Hot Intergalactic Medium, Gen & Ostriker 1999, Dave et 
al 2001). It is also possible to derive this result heuris- 
tically as follows. Let A be a wavelength of bulk mo- 
tion of the intergalactic medium in the near Universe; 
when these waves collide and break, the thermal veloc- 
ity of the shocked gas will typically be w ~ HqX, or 
V « 100/io.7(A/1.5 Mpc) km s~^ for A ~ the cluster size. 
Thus if gas heating takes place mainly at the 'nodes' of 
cluster-scale mass clumping, the thermal motion would 
involve a value of v that places the gas in the 10^~^ K 
range of temperatures. 

Since the arrival of its theoretical prediction, search for 
the WHIM has been an ongoing effort, with some suc- 
cess but no clinching evidence as yet. At the low end of 
the WHIM temperature scale, far UV absorption lines 
have been reported (e.g., Richter et al. 2008, Stocke et 
al. 2006, Tripp et al. 2006, Savage et al. 2005, Danforth 
and ShuU 2005). At temperatures where the bulk of the 
WHIM is expected to be, the detection of O VII and O 
VIII absorption lines in the spectrum of distant quasars 



(Nicastro 2005) is a more debatable result (Kaastra et al 
2006). Although soft X-ray and EUV emission in regions 
of galaxy concentration (Werner et al 2008, Mannucci et 
al 2007, Zappacosta et al 2002, Mittaz et al 1998) may 
also be the signature of warm filaments, the definitive 
proof of this interpretation, viz. an identifiable line at 
the appropriate redshift, is still not available. More pre- 
cisely, the clinching line signature of O VII was claimed 
(Kaastra et al 2003, Finoguenov et al 2003) and refuted 
(Lieu & Mittaz 2005, Takei 2008). 

We examine if limits can be placed on the WHIM from 
less direct measurements. The column density of free 
electrons in the intergalactic medium is proportional to 
the average WHIM density pwHiM in a given volume, and 
it is independent of its state of clumping. To see this, al- 
low the WHIM matter to reside in clumps of number 
density n, each having mass density p and radius r. By 
mass conservation inside a large spherical volume of ra- 
dius R, we have pwHiM ^ npr^. The average electron 
column along a random direction and for a given WHIM 
ionization fraction, being proportional to the product of 
the mass column pr of one clump and the number of in- 
tercepted clumps ^ nr^R, will then ^ pwHiM^, i-e. a 
constant for any sightline of length ~ R that penetrated 
our large volume. 

There are several ways of constraining this electron 
column. The Sunyaev-Zel'dovich effect (of the WHIM) 
will be at a very low level, as will reionization, with 
STqmb/Tcmb ~ o-ThneikT / {meC^) ~ 10"'' for character- 
istic WHIM parameters at the minimum (homogeneous) 
overdensity of 6 — fiwHiM/^^c ~ 0.02, viz. 

'fcT w 0.1 keV {T w 10^ K); 
< -e- 10-^0^) (^)^ cm"3; (1) 
^ « 2 Gpc, 

where in obtaining we assumed a fully ionized pure hy- 
drogen plasma. For comparison, the Sunyaev-Zel'dovich 
effect of a typical massive cluster (T = lO^K, Ue = lO^'^— 
lO-"* cm-3, I - 2 Mpc) is STcmb/Tcmb - - 10"^ 
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(sec Hcrnandcz-Monteagudo et al. 2008 for a discussion 
of tlie SZE from tlie WHIM). Anotlier limit might be af- 
forded by the realization that such a column can cause 
frequency dependent delay on the timescale of minutes 
to hours in the arrival of ~ 100 MHz emission from dis- 
tant quasars, except quasars with so rapid a variability 
as to avail themselves for this test arc those very ones 
affected by plasma scintillations in our local interstellar 
medium (Dennett- Thorpe and de Bruyn 2002). Angular 
broadening of quasars caused by WHIM-like scintillation 
was assessed by Lazio et al 2008, who concluded that the 
only way of securing a useful observational limit is if an 
AGN is found to 'twinkle' at a position close to that of a 
pulsar, as data about the latter will enable us to take out 
the interstellar effects of our Galaxy. Thus, the situation 
regarding these 'tangential' probes is that they too do 
not deliver any useful verdict. 

Here we take a step backwards by returning to the 
prospect of direct WHIM filament detection at the out- 
skirts of clusters, for reasons that would soon become 
clear. The most powerful argument for excluding any 
warm gas association with the central soft X-ray excess 
of clusters (e.g. Fabian 1996) is the large radiative cool- 
ing rate, since the gas (which at 10^ K is already at the 
peak temperature of its cooling curve) has to be clumped 
to co-exist with the hot cluster medium, making it radi- 
ate even faster. On the other hand, the scenario of soft 
photons seen in projection from a line-of-sight WHIM 
filament at the cluster's outskirts is more attractive, not 
only because the lack of physical contact between the 
two phases now alleviates the warm gas from its former 
problems, but also because numerical codes of structure 
formation (Gen & Ostrikcr 1999, Dave 2001, Cheng et 
al 2005) do expect such filaments to preferentially con- 
verge at clusters and groups, which are the 'knots' of the 
WHIM network. 

The purpose of this paper is to calculate the first an- 
alytical formula for the radial profile of soft X-ray sur- 
face brightness by employing a simple filament model. 
We shall see that this already affords a means of plac- 
ing rather useful observational constraints on the WHIM 
parameters at the outskirt of a cluster, including and es- 
pecially the total mass of the warm baryons. We will 
then propose the possibility of a cluster's soft excess 
being due in part to the strong lensing of background 
WHIM clumps, the largest population of which is galaxy 
groups. This is quite a radical approach to the problem: 
it may explain the lack of identifiable O VII absorption 
at the redshift of the candidate cluster, because if dis- 
tant WHIM emissions are superposed and focused by the 
cluster then their O VII absorptions would be likewise, 
leading to a smearing of the line, or even a 're-location' 
of the line to an altogether different redshift. 

2. DIFFICULTIES WITH A STRAIGHTFORWARD WHIM 
INTERPRETATION OF THE CLUSTER SOFT X-RAY 

EXCESS 

We examine more carefully why the idea mentioned in 
the previous section of the central soft excess being emit- 
ted by warm intracluster gas unrelated to the WHIM was 
not well received. Since the center of a cluster is perme- 
ated by the hot X-ray gas at (or close to) the virial tem- 
perature, the warm component can only exist for any 
conceivable length of time if it is clumped into dense 



clouds to ensure pressure equipartition between the two 
phases. Now the typical parameters of the hot gas are 
kT « a few keV and Ug ~ 10"'^ cm~'^, and because 
the pressure requirement implies equality of the prod- 
ucts TigT, we have rig > 0.01 cm~^ for these clouds. The 
immediate problem is the radiative cooling time, 

= 6 X 10« (^) ' {^;^y' years, (2) 

which means the cloiid cools and collapses on a timescale 
far shorter than the age of a cluster; it is thermally unsta- 
ble (note that in eq. (2) we omitted the contribution from 
line emission which will reduce Twarm even more). Apart 
from cooling, the cloud has difficulty in sustaining itself 
against photo-ionization by the hot gas, which occurs in 
a time Tphoto « 2 x 10''/(F/10^ ph cm^^ s~^) years for 
O VII, where F is the X-ray flux from the gas. More- 
over, any 'balance' that may result from the two oppos- 
ing mechanisms of cooling and ionization is bound to be 
extremely precarious. 

In Section 1 we also presented the cluster's soft X- 
ray excess emission as possible evidence for the WHIM. 
Among the possible models for the soft excess (e.g., see 
review by Durret et al. 2008) is the Cheng et al. (2005) 
proposal of very dense gas in pressure equilibrium as- 
sociated with merging structures within clusters. This 
model relieves the cooling problem described above in 
this section, since its non-equilibrium emission can be 
sustained for periods longer than those calculated in 
Equation 2. Alternatively, the WHIM may be seen in 
projection against the cluster's X-ray emission, although 
the filaments' brightness is approximately one order of 
magnitude too faint to explain the soft excess flux (e.g., 
Dolag et al. 2006; Mittaz et al. 2004). To date, the soft 
excess has been detected in several clusters at low red- 
shift; the Bonamente et al. (2002) blind survey of the ex- 
cess emission found that ~ 50% of the clusters have sta- 
tistically signiflcant evidence for the excess. Since clus- 
ters with high S/N observations preferentially have soft 
excess, the previous statistic is probably a lower limit. 
The excess flux is typically 10 — 20% relative to the hot 
IGM contribution in the soft X-ray band, with a trend 
of increasing relative flux at larger radii (e.g.. Lieu et al. 
1999). 

Let us now elaborate on the filament origin of the soft 
excess, by adopting as representative WHIM parameters 
those of the thermal model of Coma cluster's soft ex- 
cess, as listed under 'warm component 1' of Lieu et al 
1996, Table 3 (see also Tables 2 and 3 of Bonamente et 
al 2003) ^. We select two regions, the 0-3 arcmin and 9 
- 12 arcmin annuli, to investigate if the radial profile of 
soft X-ray surface brightness can accommodate the fila- 
ment model. For these regions the temperature is still 
kT 0.1 keV, but it is the emission integral of the opti- 
cally thin gas filament, with the value 

m^EM-A^ nlL*A = 4.79 x lO^'""^ cm"^ (3) 

for - 3 arcmin, where L* and A are respectively the 
effective length and cross-sectional area of the emitting 

^ According to Lieu et al. (1996), the warm component in the 
0-3 and 9-12 arcmin regions has a spectral normalization constant 
of 3.7 ± 1.0 X lO^'"* and 1.2 ± 0.6 X 10^''', respectively; this normal- 
ization is proportional to the EM via EM = 2.3 x 10^* x Norm x 
(area of annulus in arc minutes square) . 
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column, and EM is the emission measure (or column den- 
sity of emission), that has the physical significance. For 
a redshift of z = 0.0232 and a Hubble constant of ft. = 
0.7, we have A = 2.15 x 10^^ cm^. Hence the effective 
length of the filaments (i.e. the total length of sections of 
emitting plasma (that may belong to different filaments) 
intercepted by the line of sight) is 



L* = n-^ ■ EM ; 



10- 



Mpc. 



(4) 



Moving outwards to the 9-12 arcmin annulus, L* is 
found to fall by 3.5 times. The observed EI rises by 
2.3 times from the value in Equation [31 viz. EI = 1.10 
X 10®^ cm^^, causing L* to decrease by ~ three times 
because the projected area A of this annulus is seven 
times larger than that of the 0-3 arcmin circle. Beware 
that if the electron density of rie — 10'^ cm^"^, which 
corresponds to a rather large WHIM overdensity of (5 w 
200, is lowered, L* will easily become ^ 1 Mpc, i.e. the 
filaments will have a dimension far exceeding the cluster 
scale, placing them at 'intercluster' venues which further 
makes it harder to justify the use of (5 « 200. 

3. ANALYTIC RADIAL PROFILE OF THE SOFT X-RAY 
SURFACE BRIGHTNESS OF WHIM FILAMENTS 

The filament model, viz. linear emitters protruding 
outwards from beyond some cluster radius, can account 
for the radially rising trend of the soft X-ray excess re- 
ported by observations (Durret et al. 2008; Bonamente 
et al 2002, 2001; Lieu et al 1999). Yet the surface bright- 
ness profile of such a configuration was never evaluated 
even under the simplest scenario to see if this 'promise' 
is deliverable. We attempt to do so here, by invoking a 
model of N filaments, each of length L, cross sectional 
area A, and electron density rig, which converge to cover 
a fraction /o (by area) of the spherical surface at some 
radius R - the radius of filamentary 'footpoints' (see Fig- 
ure 1). The WHIM filling factor at any radius r is then 
given by 




/o (f ) 
for r < R. 



ioT r > R 



(5) 



The soft X-ray surface brightness of the WHIM, for a 
line of sight with impact parameter b, is then given by 
Sx (b) = EM • Aee/47r, in which Age is the WHIM plasma 
emissivity, and the emission measure 

f°° dr 
f{r)dz = 2fonlR' / 



EM 



nlL* 



-62' 
(6) 

where f — R ior b < R and f = 6 for & > _R. 

The integral of eq. (6) may readily be evaluated, yield- 
ing 



EM = <^ , 



tan" 



for b < R: 



for b> R. 



(7) 



The two results match each other smoothly at the 'foot- 
point' radius b = R {as they should), and this is also 
the radius of maximum brightness. Beyond b = R, EM 
falls off with b as EM ^ 1/6 ~ 1/9 where 9 is the an- 
gular distance away from the cluster center. Beneath 
b — R, however, EM is almost a constant, changing (de- 
creasing) only slightly from EM = TrfQ-n^R a± b — R to 




area A 



Fig. 1. — A simple model of WHIM filaments converging to some 
critical radius ij of a cluster; at this radius the surface cover factor 
by the filaments reaches the maximum value of /o mentioned in 
the text. Each filament has length L, cross sectional area A, and 
is optically thin to soft X-rays. 

EM = 2/origi? at 6 = 0. Since we understood from sec- 
tion 2 why the footpoints have to be at i? > a few Mpc 
to ensure 'segregation' between the warm and hot gases, 
and since all the observations of the rising radial trend of 
soft excesses (see references above) have thus far involved 
physical radii less than a few Mpc, this means it is the 
b < R behavior of EM that matters. The encouraging 
news is that here, the near flatness of the profile over the 
entire scale height of the virialized X-ray cluster emis- 
sion can account for the observed trend of the relative 
soft excess. 

In order to compare this model with the Bonamente et. 
al (2003) measurements of the soft X-ray surface bright- 
ness in the neighborhood of the Coma cluster, we plot 
in Figure 2 the surface brightness Sx for the parameter 
values R = 1 Mpc, Ue — 10~^ cm~'^, /o = 0.5, and with 
the emissivity being that of a plasma at kT = 0.1 keV, 
metal abundance A = 0, and averaged over the ROSAT 
1/4 keV band using the APEC model of Smith et al. 
(2001) at Aee = 4.5 x lO-^^ counts cm^ s^^ The red 
Hue in Figure 2 represents the 1/4 keV ROSAT All-Sky 
Survey background within 2-5 degrees of Coma. Given 
that the excess emission at the outskirts of Coma is on 
par with the local 1/4 keV background, we see that the 
filament model can account for the observed magnitude 
of the outer soft excess. 

Such a simple model yields no central peaking of the 
absolute surface brightness, which could be a conse- 
quence of radial dependence in rie. A, or L; although 
in this paper we would like to discuss an altogether dif- 
ferent mechanism that might also play an important role. 
Note that the effective emission length is L* ~ /qR irre- 
spective of the actual filament length L. The parameter 
L does determine the total mass of the WHIM associated 
with the cluster, which is 



MwHIM - 



47rr'^nempf{r)dr = 4.25 x 10^ 



R 



3 (f^ 

10-3cm-3y V0.5 



5 Mpc 



EM 



3 X 1018 cm 



Mr. 



(8) 



Eq. (8) constrains the filament model: there is far less 
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Fig. 2. — Intensity of the filament emission following the model 
of Section 3. The surface brightness has been multiplied by the 
average effective area A^f j ~ 150 cm-^ of the PSPC instrument 
in the 1/4 keV band (R2 band, Snowden et al. 1998), in order to 
compare this detector-dependent intensity to the value measured 
by Bonamente et al. (2003) in the neighborhood of the Coma 
cluster (shown as the red dashed line). 

room for manipulation of the parameters rie and L than 
previously thought (EM is fixed by observations), and 
arbitrary choice of their values can easily lead to an un- 
acceptable proliferation of the mass budget. 

4. THE MASS PROFILE OF A CLUSTER OF GALAXIES 

Before discussing the way in which Icnsing by a clus- 
ter's gravitational field can affect its soft excess emission, 
we first summarize our understanding of the mass distri- 
bution of clusters, which is found to follow the general- 
ized NFW profile (see Mahhavi et al 2007 for a review). 
Although in the 'authentic' version of it (Navarro et al 
1995, 1996, 1997) the density p(r) - 1/r for smaU r, 
steepening gradually to 1 /r^ at intermediate r and finally 
reaching at large r, the points at which these tran- 
sitions occur are not sharp (i.e. the index of r changes 
gradually and continuously), and are determined by the 
core radius and concentration parameter, both of which 
can vary from cluster to cluster. 

Especially for the inner radii, there is considerable de- 
bate about the form of p{r). For example, Tonini et al. 
(2006) and Schmidt and Allen (2008) find fiatter than a 
dependence, Navarro et al. (2004) and Diemand et 
al. (2004, 2005) find inner slopes consistent with the 1/r 
slope (r"" with a ~ 1.1±0.4), while Moore et al. (1998) 
determined an inner slope of a > 1.4. In any case, even 
those who favor a more gentle slope of |a| < 1 invoked 
high concentration parameter that indicates considerable 
mass inside these relatively smaller radii. Moreover, lit- 
tle is known about the nuclear region of r < 10 kpc, 
where the frequent presence of a bright central galaxy 
could cause a re-steepening of the slope. 

We shall adopt the best-fit NFW profile of the com- 
bined Subaru and Hubble ACS data of the cluster 
Abell 1689 (Broadhurst et al 2005, Diego et al 2005), 
which represent a comprehensive, model independent 
weak and strong lensing mass survey of the cluster, as 
a sufficiently representative mass profile for the purpose 
of this paper. From Table 3 of Broadhurst et al 2005 
one finds the SIS density scaling of p{r) ~ 1 /r^ between 
approximately r = rnim = 27//io.7 kpc and r — rmax = 
43O//10.7 kpc. At the r = (14 - 27)//io.7 kpc radn the 



slope has still not reached \a\ = 1, and further inwards 
are terra incognita: as already mentioned the value of 
\a\ here could increase again. The normalization for p(r) 
within the SIS zone corresponds to an enclosed mass per 
unit radius of 

= 3 X 10^5 M© Mpc-\ for n„i„ < r < w, (9) 



irrespective of the value of the Hubble constant h. 

The choice of Abell 1689 as a representative case is 
motivated by the following considerations. Abell 1689 
is a very relaxed massive cluster, according to the mor- 
phology of its X-ray emission (e.g., Xue and Wu 2008), 
with no evidence of significant substructure; moreover, 
the mass profile based on the lensig data of Broadhurst 
et al. (2005) is smooth in the region of interest, lending 
further support to its relaxed nature; finally, the lensing 
data for this cluster is of the highest quality available. 

5. GRAVITATIONAL LENSING OF SOFT X-RAY EMISSION 
FROM BACKGROUND WHIM CLOUDS 

In setting up the prerequisites let b be the impact pa- 
rameter of a light ray under the influence of some spher- 
ically symmetric cluster lens, and bo that of the same 
ray in the absence of the lens. Since the lens does not 
expand with the Hubble flow, b and bo are physical di- 
mensions independent of the epoch of lensing. Provided 
b falls within the SIS part of the cluster mass profile, i.e. 
b satisfies r,„in < b < Tmax where the two limits were 
defined in section 4, the deflection angle ip will be given 

by 



4GM 



arccos 
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(10) 

(see e.g. Lieu & Mittaz 2005). Essentially ip is a constant 
within this SIS zone: it ranges from tfj — 2ttGM/ (c^rmax) 
at r < Tmax to V = 4G'M/(c2rmax) at r < r^ax- 
Adopting the former, and applying eq. (9), one obtains 
ip = 9.3 X 10"'' which we shall henceforth approximate 
as 



= 1^0 « 10 



-3 



The lens equation 



then simplifles to 



where we defined 



b-bo 
D 



bo^b - Dtpo, 



D 



DiDu 



(11) 

(12) 

(13) 
(14) 



and the reader is referred to Figure 3 for the meaning 
of the various comoving distance indicators on the right 
side of eq. (13). 

Now consider a small element of the background source 
that maps to an element on the lensing plane defined 
by radial and transverse intervals db and bdip, with an 
area of dA — bdbdip. When the lens is removed, the 
rays from this same source element would pass through 
another element with area dAo = bo dbo dip. The area (or 
brightness) magnification of this element due to lensing 
is therefore 

^ A h 

(15) 



dA 




bdb 


dAo 




bo dbo 
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Fig. 3. — Strong gravitational lensing of a background source 
AB by a singular isothermal sphere (SIS) which deflects light by 
a constant angle i/'o- The two rays from point A that reach the 
observer after they are bent by the lens (not drawn to scale) have 
the impact parameters b = ±60 1 ^nd cast images of A that form 
part of an axially symmetric ring - the Einstein ring of physical 
radius = |bo| on the lensing plane. The impact parameter of 
the uiileiised ray that connects point B to the observer O also 
happens to be 6 = 60. After lensing, the actual rays that allow O 
to see B have impact parameters of opposite signs, viz. b = bi and 
b = b2, where |6i| = 2|6o| = 2rB. If the arrow AB were to rotate 
azimuthally about the optically axis AO to form a circular source 
disk of unlensed radius \bo\ = te on the lensing plane, the resulting 
image will also be a circular disk, but of radius 2rE on this same 
plane. Hence the brightness magnification m of the source will be 
four. If AB has smaller size than currently, then m > 4. If AB 
is larger, we will have m < 4 and the source will be too big to fit 
within the strong lensing limit. 

The modulus signs are needed because the ratio can be 
negative if caustics are involved, leading to image rever- 
sals. 

We proceed to identify the region of strong lensing, 
which has the characteristic scale length Tc given by the 
circular caustic that are images of the central spot, where 
the magnification is infinite. This limiting radius, known 

as the Einstein radius ve. is fixed by the condition 60 = 
which becomes, after applying eq. (12), 



(16) 



For finite b whilst maintaining h ~ Ve, we have from eq. 
(12) 

b^^b-Di^^^b-TE. (17) 

For b > Te the right side of eq. (15) is positive. For 
b < te, however, it is negative, and the source pixel is 
on the other side of the optical center as the image pixel, 
with values of differing by tt; images here are reversed. 
In this case we can still write 



l&ol =rE - b = rE - \b\. 



(18) 



Thus for every source pixel with |6o| < there are two 
image circles, with radii |6| = ± |6o|. As shown in 
Figure 3, a small and centrally aligned source disk of 
radius |&o| = te is imaged twice, once onto the annulus 
between |6| = and |6| — 2rE and once, reversed, onto 
the disk |6| < rE- For this source, it's clear that \db\ = 
\dbo\, so the magnifications of the two images are, from 
eq. (14), 

b 



m± = 



bo 



± 1. 



(19) 



The total magnification is m = m(&o) = m+(6o) + 
TO_(6o) = 2r£;/|6o|- Hence a centrally aligned source 
of projected physical radius |6o| = "^e, or area Aq = nr^ 
on the lensing plane and without the lens, will become 



four times brighter with the help of the lens, because it 
will then occupy the new area of 



/' 

Jo 



2Trm{bo)bodbo = 47rr| 



(20) 



Further, according to eq. (18) a smaller aligned source 
will be magnified by even more, while a much larger 
source will not be affected very much at all, i.e. it will 
only be weakly Icnscd. 

There are two limiting scenarios of relevance to WHIM 
emissions. If a WHIM source lies closely behind the clus- 
ter because it is a filament section visible to the observer, 
we will have Di k, Dg, rE ~ Digipo, and the physical size 
of the source will also be w &o- Hence the strong lens- 
ing criterion of 60 ^ te would imply a filament width 
w < DigipQ. Since Du < L where L w 5 Mpc is the scale 
length of a filament, we obtain after applying eq. (11) 
the requirement of w < 5 kpc. It would seem unlikely 
that WHIM filaments are as thin as this. Perhaps an 
even more serious problem is the smallncss of the impact 
parameters b for rays from such a source, as 5 ~ r^; « 
5 kpc also, i.e. now the lensing takes place inside the 
nuclear part of the cluster where, as already discussed 
in section 4, the density profile is not well known and 
is probably p{r) ~ 1/r, which gives rise to less lensing 
than that of the SIS profile currently being considered. 
In any case, the effect of this kind of brightening process 
is spatially very limited, it will at best enhance the soft 
excess in the central few kpc of the cluster which is not 
easily resolved by X-ray observatories. 

The more interesting second possibility is one that in- 
volves distant background WIIIM clouds, which could 
find themselves lying directly behind our optical center 
because as will be shown below the scenario under which 
any random sightline may intercept one of these clouds 
with appreciable chance is not so demanding. Assum- 
ing (reasonably) that the WHIM network extends to a 
distance ^ ~ 1 Gpc behind the cluster and the observer 
is much closer to the cluster than that, we now have 
Dis Ds '> Di, rE < Disipoi but the strong lensing 
criterion of 6o ^ rs still constrains the physical radius of 
the source to a < Distpo- For a 'mid-way' source oi Dis ~ 
O.S/i-Q^y Gpc this yields a < O.bh^j Mpc with the help of 
eq. (11). The lensing is then 'self consistent' because it 
falls within the SIS zone of the cluster, and moreover soft 
X-ray emitting WHIM clouds with scale sizes 0.5 Mpc or 
smaller are found by both simulations (see Figure 7 of 
Mittaz et al 2004a) and observations (see Soltan et al 
1996, 1997 and Finoguenov et al 2007, especially the last 
paper) to be associated with galaxy groups and galaxy 
halos. 

In fact, if all of the n « 1.56 x lO"^/!^ ,. Mpc"^ pop- 
ulation of galaxy groups of the ESO survey (Ramella et 
al 2002) possesses soft X-ray halos of 0.5 Mpc radius and 



at a density of 



10 cm this would be suffi- 



cient to account for all the missing baryons of section 
1 as WHIM constituents. The halo interception proba- 
bility by a sightline through the cluster center (or any 
sightline) is 



p = -Kna 



1=12 



n 



1.56 X 10-4 Mpc 

2 



a 



0.5 Mpc J V 1 Gpc 



% 



(21) 
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for h =0.7, and is by no means negligible as already 
mentioned. The product of p and the emission measure of 
EM « n^a per halo gives the average EM for the sightline 
in the same form as it was used it sections 2 and 3. It is 



EM = 2 X 10 



17 



n 



1.56 X 10-4 Mpc" 

3 



10"'^ cm~'^ 



0.5 Mpc 



1 Gpc 



cm 



(22) 



for h — 0.7. Without strong lensing, EM is as given by 
eq. (21) and is « 2.5 times below even the brightness of 
Coma's 9-12 arcmin soft excess, section 2. It is there- 
fore barely distinguishable from the background. With 
lensing, however, EM is enhanced by a factor « 4, eq. 
(19), and can now account for 50 % of the 0-3 arcmin 
soft excess. 

Moreover, the lenscd flux is superposed over a cluster 
radius of rimago — "^te — where is defined in eqs. 
(15) and (13). For the value of Dis = O-bh^^ Mpc we 
adopted, and using Di = 100/iq 7 Mpc as Coma's dis- 
tance, we obtain nmago ~ 0.17/i(^7 Mpc assuming a flat 
Universe where Ds — Di -\- Dig , whereas the mean phys- 
ical radius of Coma's 9-12 arcmin annulus is Q.?tlh^\ 
Mpc. Thus the flux is not dispersed into a large area 
that includes the outskirts; rather, it affects a substan- 
tial and resolvable part of the cluster's core region. This 
can then represent a non-negligible fraction of the soft 
excess there. Note that the calculation here is conserva- 
tive, because each background WHIM halo was treated 
as uniform. If their emissions are centrally peaked (due 
e.g. to a radial scaling of the filaments' density, cross 
sectional area, or length, section 3), the strong lensing 
will enhance the surface brightness of these peak regions 
by more than a factor of four, eq. (18), and the resulting 
observed brightness profile of soft X-rays will likewise be 
even sharper and more interesting. Under such a sce- 
nario, lensing of background halos could account for all 
of the inner soft excess of some clusters. 

Given the degrees of freedom of our present model and 
the uncertainty in the mass profile gradient, it is not 
possible to establish with certainty what fraction of the 
excess is due to strong lensing. Our estimate of a 10% 
probability of lensing, together with the detection of soft 
excess in >50% of the clusters, suggests that additional 
mechanisms are responsible for the excess, especially at 
large radii. Thermal emission due to dense gas (Cheng 
et al. 2005) or by the filaments themselves (e.g., Mittaz 
et al. 2004; Dolag et al. 2006), or non-thermal emission 
(e.g., Sarazin and Lieu 1998; Blasi and Colafrancesco 
1999) are viable alternatives to complement the contri- 
bution due to the lensing effect. 

6. CONCLUSION 

Although there are various indirect ways of detecting 
the WHIM and confirming (or refuting) the hypothesis 
that most of the missing baryons at low rcdshift are in 
this phase, the more promising search would still appear 
to be those involving emission and absorption of soft 
X-rays in the vicinity of clusters of galaxies where the 
WHIM is expected to congregate. 

In this paper we developed a simple model of WHIM 
filaments to account for the surface brightness profile 
of the outer cluster soft X-ray excess, in particular the 



radial rise of the fractional excess relative to the nor- 
mal X-ray emission. The inner excess, however, necessi- 
tates a separate intepretation. We propose that this is 
due to the strong lensing of background soft X-ray halos 
(most likely associated with galaxy groups) by the cen- 
tral gravitational potential of clusters. When applied 
to the Coma cluster, using reasonable parameters for 
these halos, we have the prospect of explaining ^ 50 
% of the 0-3 arcmin soft excess The calculations in 
this paper invokes filaments of constant length, cross- 
sectional area and electron density, and thus they must 
be considered as a proof-of-principle for a more realistic 
model; for example, one which allows the filament prop- 
erties to change with cluster radius. In the case of some 
clusters, a possible additional contribution to this excess 
is the strong lensing by the cluster of WHIM emissions 
associated with background galaxy groups. If a central 
gradient is present in the the emission of some of these 
aligned sources, the lensing will be particularly effective, 
and will further modify the soft excess profile of the can- 
didate cluster. 
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